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Abstract 

We present a status report of our systematic theoretical and phenomeno- 
logical study of QCD-instanton induced processes in deep-inelastic scattering. 
We show that this regime plays a distinguished role for studying manifesta- 
tions of QCD-instantons, since the typical hard momentum scale Q provides a 
dynamical infrared cutoff for the instanton size p<0(l/Q). For deep-inelastic 
scattering at HERA, we present a preliminary theoretical estimate of the total 
instanton-induced cross-section (subject to appropriate kinematical cuts). It is 
surprisingly large, in the 0(1 — 100) pb range, albeit still uncertain. We report 
on our investigation of the discovery potential for instanton-induced events 
at HERA by means of a Monte Carlo event generator. It is based on a de- 
tailed study of the characteristic signatures of the final state, like a large total 
transverse energy, Et = O(20) GeV, a large multiplicity, n = 0(25), and a 
flavour-democratic production of hadrons. A combination of event shape infor- 
mation with searches of K° mesons, muons, and multiplicity cuts might help to 
discriminate further the QCD-instanton induced processes from the standard 
perturbative QCD background. 



*To be published in: Quarks '96, Proc. IXth International Seminar, Yaroslavl, Russia, May 5-11, 
1996. 



1 Introduction 



The Standard Model of strong (QCD) and electro- weak (QFD) interactions is remark- 
ably successful. Its perturbative formulation appears to be theoretically consistent 
and agrees with present experiments. Yet, the existence of Adler-Bell-Jackiw anoma- 
lies implies that there are also processes that cannot be described by conventional 
perturbation theory ||. They give rise to a violation of certain fermionic quantum 
numbers, notably chirality (Qs) in (massless) QCD and B + L in QFD. 

Such anomalous processes are induced by instantons || which represent tunnelling 
processes in Yang-Mills gauge theories, associated with the highly degenerate vacuum 
structure ||]. 

An experimental discovery of such a novel, non-perturbative manifestation of non- 
abelian gauge theories would clearly be of basic significance. 

A number of results has revived the interest in instanton-induced processes during 
recent years: 

• First of all, it was shown || |(| that the generic exponential suppression of these 
tunnelling rates, cx exp(— An/a), may be overcome at high energies, mainly due 
to multi-gauge boson emission in addition to the minimally required fermionic 
final state. 

• A pioneering and encouraging theoretical estimate of the size of the instanton (/) 
induced contribution to the gluon structure functions in deep-inelastic scattering 
(DIS) was recently performed in Ref. 0. It was argued that it is possible to 
isolate a well-defined and sizable instanton contribution in the regime of small 
QCD-gauge coupling, on account of the (large) photon virtuality Q 2 . While the 
instanton-induced contribution to the gluon structure functions turned out to 
be small at larger values of the Bjorken variable x, it was found in Ref. [0 to 
increase dramatically towards smaller x. 

• Last not least, a systematic phenomeno logical and theoretical study is under 
way || D, [10], |TJ, [TJ] , which clearly indicates that deep-inelastic ep scattering at 



HERA now offers a unique window to experimentally detect QCD-instanton in- 
duced processes through their characteristic final-state signature. The searches 
for instanton-induced events have just started at HERA and a first upper limit of 
0.9 nb at 95% confidence level for the cross-section of QCD-instanton induced 
events has been placed by the HI Collaboration fl3f . New, improved search 
strategies are being developped [ I 1] with the help of a Monte Carlo generator 
(QCDINS 1.3) |§ for instanton-induced events. 

In this review, we report on the present status of our broad and systematic study of 
QCD-instanton induced processes in deep- inelastic scattering || £| [II], [IT], [T2], |T5j]. 



In Sect. |2|, we start by summarizing the important theoretical result [E2] that 
instanton-induced processes in deep-inelastic scattering do not suffer from the usu- 
ally encountered infrared (IR) divergencies associated with the integration over the 
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instanton size. In fact, the inverse hard momentum scale in DIS provides a dynamical 
IR cutoff for the instanton size. We point out the close analogy of QCD-instanton 
induced processes in DIS to those in QFD. In Sect. ||, we present the current status 
of our ongoing theoretical study |14| of the /-induced contribution to deep-inelastic 
lepto-production. Of interest is, of course, the first (preliminary) estimate of the 
total /-induced cross-section at HERA. Section f| is devoted to our phenomenologi- 
cal investigation of the discovery potential for instanton-induced events at HERA by 
means of a Monte Carlo event generator |5], The emphasis rests on the study of 
the characteristic signatures of the final state. We also report on the first searches 
of instanton-induced events in DIS at HERA. Finally, we discuss improved search 
strategies, which can serve to enhance the signal from /-induced events relative to 
the background from normal DIS events. 



2 What is Special in Deep-Inelastic Scattering? 



As is well known, instanton calculations in QCD are generically plagued by IR diver- 
gencies associated with the integration over the instanton size ||16|| . In this Section, we 
shall demonstrate explicitly that the typical hard momentum scale Q in deep-inelastic 
scattering provides a dynamical IR cutoff for the instanton size p I <0{1/ Q), such 
that the integration over the size parameter is finite |12[]. Let us consider the ampli- 
tude for the relevant instanton-induced chirality-violating photon-gluon process (see 
Fig.0), 



i* + § =>■ \& + <ih] + n § s ; 

flavours 



[AQ 5 = A(Q R -Q L ) =2 n/) 
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The strategy is to first set up the respective Green's function according to standard 
instanton-perturbation theory in Euclidean space [EL O, [TR [Tjl FJ, then to Fourier 
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transform to momentum space, to amputate the external legs, and finally to continue 
the result to Minkowski space. The basic building blocks are (in Euclidean space and 
in the singular gauge): 



i) The classical instanton gauge fieldQ A^/ 

2tt 2 



, , . <y x - ay \ t f 



(2) 



n, 



X 2 ' 



(3) 



depending on the various collective coordinates, the instanton size pi and the colour 
orientation matrices Uj. The Uj matrices involve both colour and spinor indices, the 
former ranging as usual only in the 2x2 upper left corner of 3 x 3 SU(3) colour 
matrices. 



ii) The quark zero modes 0, k and 0, 



1 



2tt 2 x 4 ifj 2 
1 



7 X 



2tt 2 x 4 n^ /2 



and 



(4) 
(5) 



iii) the quark propagators in the instanton background [Ij] 

S (I) (x,y) = 



U x U y 



x-y 



2ir 2 (x - y) 4 



1 + P 2 i 



Ujxyllj 



X 2y2 



(6) 



Uix (x - y) a„ y U\ 



An 2 x 2 (x — y) 2 y 4 H y 



S (I \x,y) 



x-y 



2n 2 (x - y) 4 



1 + P? 



UjxyUj 



x 2 y 2 



(7) 



p] °» 



Uix (x -y)yU\ 



4tt 2 Ii x x 4 (x - y) 2 y 2 



1 We use the standard notations, in Euclidean space: oy, = (— ia, 1), = (i<j, 1), and in 
Minkowski space: er^ = (1,<?), er^ = (1,— <y), where a are the Pauli matrices. Furthermore, for any 
four- vector w M , we use the shorthand v = v ■ a, v = v ■ a. 
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Figure 2: Instanton-induced chirality-violating process, 7* + g — > qTJ + q_R + n g g, 
for one massless flavour (rif = 1), in leading semi-classical approximation. The cor- 
responding Green's function involves the products of the appropriate classical fields 
(lines ending at blobs) as well as the quark propagator in the instanton background 
(quark line with central blob). 



For simplicity, we concentrate here on the case of one flavour (n/ = 1 in Eq. ([]]) ) , 
with the generalization to a larger number of flavours being straightforward. 

The relevant diagrams for the exclusive process of interest, Eq. ([!]), are displayed 
in Fig. in leading semi-classical approximation. The amplitude is expressed in 
terms of an integral over the collective coordinates pi and the colour orientation Uj, 



% (7 + g -> Ql + Qfl + n g g) 



[duj [^fd( Pl ^ r )A; ai - an °( Pl ,u 

J J Pi 



l (0 



1,2,3, 



where 



d(pi,fi r 



d 



2tt 



exp 



2tt 



a s (p r 



(pi 



(9) 



a s (p r 

denotes the instanton density @, fL9| , with \i r being the renormalization scale. 
The density (|), with the leading-order expression for a s (p r ), 

4 Tj- 2 

a s (p r ) = - - /<<2N ; /3 = ll--n/, (10) 



satisfies renormalization-group invariance at the 1-loop levelQ. The constant d is given 
by 

2 Two-loop renormalization-group invariance of the density (^J) can be achieved pC)| by replacing 
in Eq. (||) /3 by /3 + fj(/?i — 12/3q), with next-to-leading order expression for a s (fi r ). 
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with C x = 0.466, C 2 = 1.54, and C 3 = 0.153, in the MS-scheme. 

Before analytic continuation, the amplitude A^ entering Eq. (|8|) takes the following 
form in Euclidean space, 



.aai...a n 

sA.fi 



-ie q lim p 2 tr \a a e g (p) ■ A {I) (p)] TT lim p 2 tx \a ai e g {pi) ■ A {I \- Pl )} x (12) 



X R (k 2 



lim (i^aX-^V^^-ArO+V^Cg,-^) lim (-i 



Xl(& 



with contributions from the diagrams on the left and right in Fig. respectively, 



cr„ lim S^ 7 - 1 (x, — fcj) (— ikx) 

kj^O 



(13) 



lim (ik 2 ) S (— A; 2 , x) cr^ «(:r) 



(14) 



and generic notation for various Fourier transforms involved, 

f(...,k,...) = jr/'re /( ,:...). (15) 

The four- vector e 9/Lt / denotes the gluon polarization- vector, whereas Xl,r are Weyl- 
spinors, satisfying the Weyl-equations, 



and 



k XL {k) =0; k XR {k) = 0, 
XL(k)x\(k) = k; Xn(k) Xn(k) = k ■ 



(16) 



(17) 



The LSZ-amputation of the classical instanton gauge field A^r in Eq. (|i~2|) and the 
quark zero modes k and <p in Eqs. (p~3T) and (0), respectively, is straightforward |||12]. 
On the other hand, the LSZ-amputation of the quark propagators and in 
Eqs. (|I~3"D and ([14]), respectively, is quite non-trivial and involves a lengthy calculation. 



For the details, we refer the interested reader to Ref. ||12|| . Here, we only present the 
final result for the scattering amplitude (W) in Minkowski space, 

T a ai ...a ng ^* +g ^— +qR + n ^^ =i j J dpjdip!, fl^pT" 



X tr 



X 



a a Uj [e g {p) ■ p - e g {p) p) Uj] ]J tr [a* Uj [e g ( Pi ) pi - e g ( Pl ) ■ Pi ) Uj 

i=l 

{ [Ui xi^a) e ] [ e V (Q. k ^ Pi) w v Xi(ki) Uj 
UiXR(k2)a lx V(q,k 2 ;pi)e\ [eXL(fa) U\ } , 
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where the four-vector V\ reads 



V\{q,k;pi) 



(q- k )\ , fc A 
-(q-k) 2 2q-k 

kx 



pJ-iq-kfK, pJ-(q-k) 2 ) (19) 



2q ■ k 

At this stage of our instanton calculation, the constraint arising from electromag- 
netic gauge-invariance, 

g"r; ai " ,s =0, (20) 



can easily be checked. Upon contracting ( |TBD with g M , it is easy to verify Eq. 
with the help of the Weyl-equations (|16D, the relations 

VpPv + VvVfj, = 2g^u, (21) 

of the a-matrices§ in Minkowski-space, and the on-shell conditions k\ = k 2 = 0. 

Let us concentrate now on the main objective of this Section: The integration 
over the instanton size pi in Eq. (|I5D is finite ||12j| . In particular, the good infrared 
behaviour (large pi) of the integrand is due to the exponential decrease of the Bessel-K 
function for large pj in Eq. (|19|), 

MQpt) Qp ^°° sj\^ exp [-Qpx] • (22) 

Thus, the typical hard scale in deep-inelastic scattering, 

Q = min {Q = f^f, yl-tq-ktf, y/-(q - k 2 )^ (> 0) , (23) 

provides a dynamical IR cutoff for the instanton size, pi<0(l/ Q) (at least in leading 
semi-classical approximation). Therefore, deep-inelastic scattering may be viewed as 
a distinguished process for studying manifestations of QCD-instantons. 

After the integration over the instanton size, which can even be performed ana- 
lytically [0, we find that the /-induced amplitude flI8| ) is well-behaved as long as 
we avoid the (collinear) singularities which arise when the internal quark virtualities, 
~ (<? — ki) 2 (c.f. Fig. H (left)), or — (q — k 2 ) 2 (c.f. Fig. || (right)), vanish. Thus, 
like in perturbative QCD, fixed-angle scattering processes at high Q 2 and moderate 
multiplicity are reliably calculable in (instanton) perturbation theory ||12|| . 

Note that the hard scale Q plays a very similar role as the vacuum expectation 
value of the Higgs field in electro- weak instanton-induced (B + L)-violation 0, ||, ||. 
Another close analogy refers to the fact, that, similar to electro- weak instanton- 
induced (B + L)-violation || ||, also QCD-instanton induced processes in DIS are 
dominated by the multiple production of vector bosons (gluons in the case of QCD; 
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W's and Z's in the case of QFD). This can be seen directly from Eq. (0), since each 
additional gluon in the final state gives rise to a factor of (ix z / o> s ) l l 2 in the amplitude. 

In analogy to electro- weak (B + L)-violation |J, one expects [0, || the sum of 
the final-state gluon contributions to exponentiate, such that the total /-induced 7*g 
cross-section takes the formQ (at large Q 2 ), 



a 



' " f x, Q 2 ) = a rgn g (x, Q 2 ) ~ / dx' I . . . -J- exp 



1 j , V dQ' 2 1 



47T 



a s (Q>) 



,(24) 



where the so-called "holy-grail function" Fix') (normalized to F(l)=l) is expected 
to decrease towards smaller x f , which implies a dramatic growth of cr^) g {x, Q 2 ) for 
decreasing x. 

In principle we could verify (^) by exploiting the exclusive amplitudes (|1|) fur- 
ther, i.e. by performing the collective-coordinate integrations, taking the modulus 
squared, integrating over the appropriate phase spacef^] and finally summing over n g . 
We will, however, adopt a much more powerful method in the next Section which, 
moreover, will allow us to go beyond the leading semi-classical approximation about 
the instanton (c.f. Sect. 



3 Structure Functions and HERA Cross-Section 

In this Section we want to concentrate on a first, preliminary estimate of the rate for I- 
induced events at HERA |14j (subject to appropriate kinematical cuts). Let us recall 
that, in the one photon exchange approximation, the unpolarized inclusive lepto- 
production cross-section is parametrized in terms of the familiar structure functions 
F 2 and F L , 



d 2 



a 



Ana 2 



yB 3 + y -f)F 2 (x B3 ,Q 2 )-^fF L (x B3 ,Q 2 



(25) 



where is the center-of-mass (cm. 
Bjorken variables are defined by 



energy of the lepton-hadron system. The 



IP -q 



2/Bj 



P ■ q 
P-k' 



(26) 



where P (k) is the four-momentum of the incoming proton (lepton). 

3 The lower IR cutoff /iy in the Q' 2 integration in Eq. (|24|) serves to regulate the (collinear) 
divergence mentioned above and plays the role of a factorization scale. 

4 This programme has been performed for the simplest /-induced process, corresponding to 
71/ = 1, n g = 0, in Ref. (La. The straightforward generalization is in progress Oj. 
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Figure 3: Imaginary part of the //-contribution to the forward 7*g scattering am- 
plitude. 



The /-induced contribution to the nucleon structure functions, e.g. F 2 (xej, Q 2 ), 
is described J7|, [8|, [14| in terms of the standard convolution of parton structure func- 
tions, e.g. ^,(2, Q 2 ), with corresponding parton densities, e.g. / s 



(27) 



Let us concentrate next on the calculation of the (dominating) /-induced contri- 
bution to the gluon structure functions |T3 . 



3.1 /-Contribution to the Gluon Structure Function 

The calculation rests J7|, |14| on the optical theorem combined 5 with the //-valley 



method |22|. [23|. 23]. The optical theorem for the virtual 7*^ — > 7*5 forward amplitude 
represents a convenient method to perform implicitly the summation over the I- 
induced multi-particle final state. 

Accordingly, the /-induced contribution to the gluon structure functions can be 
obtained by simple projections from the imaginary part of the //-contribution to the 
forward 7* g scattering amplitude (see Fig. 



W« (p,9) 



7T 



lmT g %\p,q) 



(2? 



-g^ +Qx 



(29) 



. p-q. 

5 For a combination of the optical theorem with the //-valley method in the context of electro- 
weak instanton-induced (B + L)-violation see Ref. pl[. 
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Figure 4: Structure of the imaginary part of the //-contribution to the forward 7*g 
scattering amplitude, in leading semi-classical approximation about the //-valley and 
at large Q 2 . 



p-q 



(30) 



where 



x = 



2p ■ q 



(31) 



is the Bjorken variable of the 7* g subprocess. 

In order to extract the required //-contribution (c.f. Fig. [5]) from the for- 
ward 7*g scattering amplitude, one first evaluates the path-integral expression 
for the corresponding Green's function in Euclidean space by expanding about 
the //-configuration (II "valley"), which is defined via the so-called "valley 

method 6 " [22, 23, 24]. After Fourier transforming to momentum space, one has to 



amputate the external legs and to analytically continue to Minkowski space. Finally, 
the imaginary part is taken. 

Along these lines, it is possible to show [EH] that, in leading semi-classical approx- 
imation about the //-valley and at large Q 2 , the dominant projection contributing to 
Tig{%i Q 2 ) has the following structure (see Fig. |]), 



g)~E <~ [ dW +) ((? - q'f 



Q 2 (p-q 



A2 



{{p + q' 



,2^3/2 



'{¥. 



JduJ^j-J^fJ *RM{p hPT ,&V) exp 



dpj 

4 



47T 



(32) 



For any fixed values of the II collective coordinates n,r = {pi, pj, U}, the pair configuration 



A^ 1 ^ is required to minimize the action within the subspace orthogonal to dA 1 ^ 1 ' jdT{. 



(ii) 
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x exp [i (p + q') ■ R\ K x (Q' Pl ) K x (Q' p~) {p lPl f 2 I — 



l 6 \ n f-2 



where Q = —q . The second projection PO), entering J r Lg (x, Q ), may be evaluated 



analogously Jl4|]. It turns out to vanish in the Bjorken limit. 

In analogy to the single instanton case discussed in Sect. ^ the amplitude (^) is 
expressed in terms of an integral over the appropriate II collective coordinates, the 
/(J)-size parameters, Pi(j), the distance vector between I and J, and the relative 
colour orientation, U. Due to conformal invariance of the classical Yang-Mills action, 
the //-valley action S^ n \^,U) depends on the distance and the /- and /-sizes 
only in a conformally invariant manner 



-R 2 + ie R + pj + pf 
PiPj 



(in Minkowski space) 



(33) 



The //-valley action itself, for the most attractive relative orientation, Uq 



y/-R 2 +ieR Q 



is given by |^3 



S {ll) (Z,Uo). 



(34) 
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96 



+ 



48 



no my no 



[3/(0 + 8] In 



2f(/(0+4) 



e+ 



^ _ 4 ^ _ 4 _ 



(35) 



In Euclidean space, where —R + ie Rq in (|33| ) is replaced by R , the //-valley action 
£^(0 corresponds to a smooth interpolation between a widely separated, weakly 
interacting //-pair configuration (S^ n \^ — > oo) = 1) and a strongly overlapping one, 
annihilating to the perturbative vacuum (S^ 77 **^ — > 2) =0). 
Finally, the combination 



M(pi,pj,£,U) exp 



47T 



a s (/i r 



-S (/7) (£,t/) 



(36) 



2tt 



a s (/x r 



12 



exp 



47T 



a s (/x r 



appearing in the //-induced amplitude (|32"D, denotes the (one- loop renormalization- 
group invariant) //-density^ 

After having defined the quantities entering Eq. (|3~2"1), let us come now to its 



physical interpretation. Complementary to the result quoted in Ref. 0, our result 



7 There are corrections to Eq. (|3^) which vanish for large £. Fortunately, we need the density 
only at sufficiently large £ (see below), such that Eq. (pq) is approximately valid. 
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yields a momentum-space picture of the //-contribution to the imaginary part of the 
forward "f*g amplitude. As is clear from the discussion in Sect. |2], the £( + )-function 
as well as the Bessel-K functions appearing in Eq. ( j32"D are easily understood as 
originating from the discontinuity across the handle in the //-induced "handbag" 
diagram in Fig. |j. The factor exp [i (p + q') ■ R], on the other hand, arises from the 
external gluons g(p) and the internal, virtual quarks q*(<?') in the // background. It 
is then tempting to introduce the notion of an I-induced q*g subprocess, along with 
its associated total cross-section [lH 



a. 



19 



x',Q 



'2\ 



8 7T 5 

3 a. 



Q' 2 (p ■ q 

i(p + q' 



A 2 



\2\ 3 / 2 



(37) 



oo oo 

duf^ fig. I ,l i RMl,>j. f > T .i,.U)^ 

o 



4 



47T 



x exp [i (p + q') ■R]K 1 [ JQ* pA KA JQ" pj {p lPl Y /z - 



y~2 



Upon partially performing the integration over q'^ in Eq. (p2[), one may easily show 
that, for large Q 2 , the /-contribution to the gluon structure function Tig-, Eq. (|29|) . 



can be expressed as JT4 



rj8(x, Q 2 )~J2 e l x \%- (g^a x ~) f d Q' 2 °&V> Q' 2 ' 



(3J 



where 



Q' 



* = -r-; (39) 

2p ■ q 

denotes the Bjorken variable of the q*g subprocess. Thus, in leading-order semi- 
classical approximation about the //-valley and at large Q 2 , the /-contribution to the 
gluon structure function Ti g has the form of a convolution of a "splitting function" 
in the /-background, 

1 ' (40) 



p (i) (£ 

q*h* \ x > 



i 7T 3 X' 



with a total cross-section (|3~7D for the /-induced subprocess, containing the essential 
instanton dynamics. A generalization of this convolution structure M to less inclusive 



quantities, e.g. single-particle inclusive distributions, is in progress [14 . 



3.2 The Total Cross-Section of the /-Induced q*g Subprocess 

In this Section we shall concentrate on the most important building block entering 
Eq. (|38|), the total cross-section (|37j ) of the /-induced q*g subprocess (c.f. Fig. |j). 
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The collective coordinate integrals in Eq. fl3"T|) can be done in the saddle-point 
approximation, the small parameter being a s (p r ). The expression to be extremized 
is the exponent in Eq. (|37p. 

47r — 

Tip, q'; PI , pj, R„ U) =i(p g + q ').R- Q' ( pi + pj) —S^^, U) , (41) 

®-s \Pr ) 

where we have used the asymptotic form (|22| ) for the Bessel-K functions, anticipating 
that, for small a s (p r ), the dominant contribution to Eq. ([37]) will come from the 
region 

q' P/(7) »i. (42) 

The corresponding saddle-point equations are most easily solved in the q*g c. m. 
system, where the solution is given by p3[ p5f] , 

4 



2 + 



Ul' 



(43) 



Pi = Pj 



16 vr 1 



(44) 



R* 



\ 



e-2 



(p*i 



pipj 



(45) 



U* 



R* 



V^R^' 

Here a;' denotes a scaling variable corresponding to the q*g c. m. energy, yi' 



Q 1 



x' 



(46) 



(47) 



Note, that the saddle-point for p, Eq. (|44|), is consistent with our ansatz ([42]), as long 
as a s (p r ) is small. 

The exponent T, Eq. (0), evaluated at the saddle-point, Eqs. ([43|-[46]), takes a 
particularly simple form [23, 25], 



T&q^p^frR^U* 



47T 



a s (p r 



(4* 



It depends only on the scaling variable x' (or, equivalently, on uj'). 

Now we are ready to perform the integrations about the saddle-point. The final 
result for the cross-section (^), for small a s (p r ), reads |L4|], 



<Jq*g{x , Q 



47T 



21/2+2 A)F(x') 



Q' 2 \a s (/i r 



exp 



47T 



■jr F(x'^ 



(49) 
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Figure 5: The holy-grail function and the valley action. 



where the "holy-grail" function F(x') is here identified as (see Eq. (fl~8D and Fig. ||), 

F(x') = S {11) (C(x')) . (50) 

Moreover, we find lO!, 



E(x') = d 2 v / 6 2 n /- 4 7r 15 / 2 - 



X 



11/2 



1 -x' 



371/ +2 



1 - X') 7 \1+X' 



X{1 ~ X) — 



2[3 F(x')-5/2 

(51) 



Let us investigate our result for the holy-grail function F(x') : Eq. (|50|), for small 
g*g cm. energy y/s', i.e. small a;' (c.f. Eq. (|47|)), 



F(x') = 1 - - 



3 /l -x'\ 2 3 (I -x'\ 3 



+ r 







1 - a/' 



(52) 



The origin of the different terms on the rhs. of Eq. (|52|) is illustrated in Fig. [| The 
first term ('t Hooft tunneling factor ||) corresponds to a purely fermionic final state. 
The next term arises from the summation of the leading-order gluon emission (c.f. 
Section 0), and the third term originates J26| from the summation of interference 
terms between the leading-order gluon emission and the gluon-propagator correction, 
see Fig. ||. 

The //-valley method allows to extrapolate (|5^) smoothly beyond instanton- 
perturbation theory, via the identification of the holy-grail function with the //-valley 
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Figure 6: Low-energy expansion of the total instanton-induced subprocess cross- 
section. The three different graphs on the rhs. correspond to the respective three 
different terms in the low-energy expansion of the holy-grail function (]52f) . 



action at the saddle point, Eq. fl50P . As we have illustrated in Fig. [| it effectively 
sums up the gluonic final-state tree-graph corrections to the leading semi-classical 
result^. However, it has been argued |27j that some initial-state and initial-state - 
final-state corrections exponentiate as well and might give rise to additional correc- 
tions, which show up in the low-energy expansion of the holy-grail function at 
order O{oj' b ). Nevertheless, it appears reasonable to trust our result (f49[-[5T|) down 
to x' = 0.2, where F(0.2) = S( /7 )(f (0.2)) ~ 1/2 (see Fig. |), a value sometimes 



advocated |28[ as the lower bound for the holy-grail function. 

As it stands, our result ( |49"D for the /-induced total q*g cross-section strongly 
depends on the renormalization scale /i r . Possible improvements, partly along the 
lines indicated in footnote ^|, will be discussed elsewhere |T^]. The main /^-dependence 
in Eq. (pi), 



4tt 



2/3 F(x') 



(53) 



has its origin in the instanton-size dependence of the //-density fl36|), (p/pj/i^) /3 ° 5(I/) , 
taken at the saddle point (f44|). A natural choice [[5], ^] for the renormalization scale 
is fi r = 1/p*. A closely related choice, which we shall adopt in the following, is to fix 
the renormalization scale, 

Mr = /i(Q') , (54) 



by the requirement, 



47T /i(Q') 



1. 



(55) 



For a formal proof, see first reference in Ref. 
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Q [GeV] Q' 2 [GeV 2 ] 



Figure 7: Left: The inverse running coupling at the renormalization scale [i(Q), 
acting as the large parameter in the saddle-point approximation. Right: The total 
cross-section of the /-induced q*g subprocess, a q * g {x' , Q /2 ), for fixed x' = 0.2. 



With this choice, the cross-section reads 0, [14]], 

\ 21/2 



a 



99 



x', Q 



4:71 



* „ ( 

Q' 2 UMQ')). 



exp 



47T 



OCsW)) 



F(x') 



(56) 



for x' > 0.2 , Q' > 5 GeV . 

The restriction on x' results mainly from our insufficient knowledge of the holy-grail 
function F(x') for small x', whereas the restriction on Q' is due to the requirement 
that the parameter An / a s (fi(Q')) should be large (see Fig. [7| (left)) in order to justify 
the saddle-point approximation. 

Note the following important feature of aqO g (x', Q' 2 ), Eq. (|56"D, as a function of Q' 
(see Fig. ^ (right)): The Q' dependences from the high inverse poweif] of a s and the 
exponential in Eq. ( |5T)D compete to produce a strong peak (maximum) far away from 
the IR region, at a new hard scale, 

A<V) = 2/3o - J" F 1 a ex P 1 1 a I A » A ' ( 57 ) 

l + /3 F(x') (l + /5 F(x')J 

much above the usual QCD scale A. This implies that the /-contribution to the gluon 
structure function J^gfe, Q 2 ), which involves an integration of (jql g (x', Q' 2 ) over Q' 2 

9 The large power 21/2 mainly reflects the number 12 of II zero modes minus 1/2 times the 
number of saddle-point integrations. 
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(c.f. Eq. (II)), 



Fg\x, Q 2 ) ~-£<%x£^ E(*' 



(5? 



Q 2 ^ 



dQ 



12 



Air 



21/2 



Q' 2 \a a (fj,(Q')) 



exp 



47T 



<KQ')) 



is asymptotically dominated by this peak and hence Q independent (scaling) in the 
Bjorken limitQ, 



The predicted approach to scaling 

Q 



/2l\ Fix') + -ft 
1 A M " (59) 



i?( x /)23/2 • 



„ d J-~2 g i x j Q 2 ) Q lar S e 



(60) 



S(x) 



tfct 



Po In 



A 2 



-2/3 F(x) 



A 2 

Q 2 



A) ^(x) 



resembles a "fractional twist" term, where the twist is sliding with x: the scaling 
violations vanish faster with increasing x. 

Note, however, that the above results for the gluon structure function depend 
to some extent on our choice of the renormalization scale fi r = fi{Q'), Eqs. (|54]), 
(|55|) . A thorough investigation of the inherent renormalization- and factorization- 
scale dependences is presently under way |TJ[. Then, we may also compare in more 
detail with the results of Ref. and sort out possible differences. 



3.3 HERA Cross-Section 

In Fig. |8] we present the resulting /-induced total ep cross-section for HERA, subject 
to the following cuts: 

• x' > x' min = (0.2, 0.3) (c.f. discussion in Sec. |37^) in the integration of Eq. (|58|); 

• XBj > a^Bj min, 2/Bj > 0.1 in the integration of Eq. (f25|). 

10 Note, that the simplest /-induced (exclusive) process, without gluons in the final state (c. f. 
Sect. ||), already gives a scaling contribution to the gluon structure function T^gix., Q 2 )- 



17 




Figure 8: Instanton-induced total ep cross-section for HERA (preliminary) with 
various cuts as indicated. 



So far, only the (dominating) gluon contribution has been taken into account. 

Apparently, the cross-section is surprisingly large. Let us, however, empha- 
size again that the inherent uncertainties associated with the renormalization- and 
factorization-scale dependences may be considerable [13]. Therefore, Fig. [5] is still to 



be considered preliminary. Furthermore, there is, of course, also a strong dependence 
on the precise value of the QCD scale A (here taken to be 234 MeV). 



4 Monte- Carlo Simulation for HERA 

A Monte-Carlo generator for QCD-instanton induced events in DIS, QCDINS 1.3 
based on HERWIG 5.8, has been developedF] [§, |i~5| . It has the following features 
built-in, which are characteristic for the underlying instanton mechanism: 

• In its cm. system, q'+p = 0, the /-induced multi-parton production is supposed 
to proceed isotropically 0, |5], [7], [TtJ. We may imagine a "fireball" in S'-wave 
configuration, decaying into gluons and at least 2nf — 1 quarks. 

• The total parton multiplicity associated with the J-subprocess is large || ||, [7], 

11 Until the final write-up in Ref. |l5| is completed, Ref. Q may be consulted as a qualitative 
guide for QCDINS. The Monte-Carlo event generator for electro- weak instanton-induced processes 
(HERB VI) in Ref. |§ shares a number of similarities with QCDINS. 
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Figure 9: Lego plot (77i a b, </>iab, -&r[GeV]) of a typical /-induced event in the HERA- 
lab system at XBj = 10~ 3 . 



(n q+g ) {1) = O(10) at HERA. (61) 



The multiplicity distribution at the parton level is taken to be a Poisson distri- 
butional ft; 



• All light flavours, including strangeness and possibly charrn^], if kinematically 
allowed, are democratically produced f2|. 

At present, we consider the characterization of /-induced events by these final-state 
features to be more robust than predictions based on cross-section estimates. 

The typical event in Fig. |9| from our Monte- Carlo generator |T5| illustrates these 



built-in characteristics: A current-quark jet along with a densely populated hadronic 
"band" of width Ai] = ±0.9 in the (r/i a b, 0i a b)-plane |§. The band reflects the isotropy 
in the /-rest system. The total Et = O(20) GeV is large as well as the multiplicity, 
^band = 0(25). Finally, there is a characteristic flavour flow: All (light) flavours are 
democratically represented in the final state. Therefore, strongly enhanced rates 
of /r°'s and /z's (from strange and charmEl decays) in the hadronic band represent 
crucial signatures for /-induced events. 

12 Strictly speaking, this holds only in the Bjorken limit O], 

13 The predictions associated with the /-induced production of charmed quarks are less certain, 
since for m q p* > 1 additional suppression factors may arise. 
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Figure 10: Numbers of K° mesons per event with 0.25 < p\ < 4.5 GeV 2 as a function 
of the HERA-lab pseudorapidity 77 for the kinematic regime 10 < Q 2 < 70 GeV 2 , 
10~ 3 < XBj < 10 -2 , 0.1 < < 0.6. Left: Predictions from various models of (non- 
diffractive) DIS (broken lines) and from QCDINS (continuous line) for /-induced 
events only. Right: HI data (points) with the fit result (dotted line) and the fraction 
of instanton induced events (solid line)), / = 0.006, obtained as described in the text. 
Both figures taken from Ref. |T3| . 



4.1 A First Search for /-Induced Events at HERA 

If a significant proportion of DIS events at HERA were induced by QCD-instantons, 
large changes in the strangeness composition of the hadronic final state would be 
expected in the "band" region, as illustrated in Fig. |iy (left). This feature of I- 
induced events has been exploited by the HI collaboration |13| who searched for an 
excess in the K° rate in the barrel of HI. 

They fitted their measured rate of K° production, shown as a function of the 
HERA-lab pseudorapidity rj in Fig. [II] (right), with various DIS models, allowing for 
a fraction / of /-induced events (see Fig. [1(] (left)). The largest value of / was obtained 
using QCDINS with LEPTO [|| and was / = 0.006 ±0.030. The values / found with 
different DIS models were all found to be consistend with being zero. Hence, a 95% 
confidence level upper limit of 0.9 nb was placed |13| on the /-induced cross-section at 



HERA in the kinematical region considered (c. f. Fig. |TU|). Note that this limit is less 
than an order of magnitude above our (preliminary) theoretical estimate presented in 
Fig. ||, which even contains stronger additional cuts (in particular x' > x' m \ n = 0.2). 
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Figure 11: Left: Et distributions in the 7*p cm. system for normal DIS (left) and 
/-induced events (right). Right: Correlation between 2E out and the i-subprocess 
c. m. energy Wj = xfs 1 . The primes indicate additional cuts in 77 to minimize next- 
to-leading-order perturbative QCD effects. Both plots refer to the kinematic region 
0.001 < x Bj < 0.01, 0.1 < y Bj < 0.6 and 20 < Q 2 < 70 GeV 2 . 



4.2 Improving the Search Strategies 



Let us finally mention some recent efforts [11] to improve the sensitivity to /-induced 



events by adding-in characteristic information on the event shapes. The first step 
consists in boosting to the ^*p cm. system and looking for events with high Et 
(c.f. Fig. |TT] (left)). We note that in this system (1+1) and (2+1) jet0 perturbative 
QCD processes deposit their energy predominantly in a plane passing through the 7*p 
direction. In contrast, the energies from /-induced events are always distributed much 
more spherically (isotropy in the /-rest system!). Therefore, one may substantially 
reduce the normal DIS background by minimizing (on an event-by-event basis) the 
quantity 

n 

E ont = nun ^2\pk-i\, (62) 

' k 

by choice of the unit vector z, normal to the 7*p direction. For standard (2+1) jet 
events from boson gluon fusion, E out is then only of order of the jet widths. In 
contrast, for /-induced events, E out ~ y/s' r /2 is large (see Fig. [□] (right)). Note, that 
the peaking of the E T and \fs' = Wj distributions in Fig. |TT] around 15 GeV directly 
reflects the shape of a q * g as function of Q' 2 oc s' (c. f. Fig. |7| (right)). The quantitative 
results from the Monte-Carlo simulation, subject to additional cuts in rj which are 



14 



As usual in DIS, the +1 refers to the proton remnant. 
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Figure 12: E ont vs. E in distributions in the 7*-proton cm. system for normal DIS 
events (top) and /-induced events (bottom) in the kinematic region 0.001 < x < 0.01, 
0.1 < y < 0.6 and 20 < Q 2 < 70 GeV 2 . The primes indicate additional cuts in r\ to 
minimize next-to-leading-order perturbative QCD effects. 



to minimize higher-order perturbative QCD effects, are displayed in Fig. |12| They 
fully confirm the qualitative expectations. The power of an /-induced event selection 
based upon event shape is certainly apparent. 

The production of large numbers of partons in the /-subprocess also leads to large 
charged particle multiplicities in the hadronic final state (c.f. Fig. |9|). This is of 
particular interest for an instanton search, if the HERA detectors can extend their 
multiplicity measurements to include the region covered by the forward trackers. 

Apart from the enhancement of the K° production rate, a second feature associ- 
ated with the flavour democracy of /-induced events is the large number of muons 
they contain. These result largely from the decay of charmed particles^. Unfortu- 
nately, their energies are rather low, with the transverse momenta of the muons in 
the laboratory frame being typically less than 1.5 GeV. Hence, their detection is a 
challenging task for the experimenters. 

Another challenging task is to study, how well the /-subprocess variables s', x' may 
be reconstructed/restricted from the hadron momenta in the "band" and the current 
jet. This would enable us to compare with the theoretical estimates of the production 
rates which, as explained in Sect. |3.2| , are only available in the range x'>0.2, Q'>5 
GeV. 
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5 Conclusions 



The experimental discovery of QCD-instanton induced events would clearly be of 
basic significance, since they correspond to a novel, non-perturbative manifestation 
of QCD. In addition, they would also provide valuable indirect information about 
(B + L)-violation in the multi-TeV region, induced by electro-weak instantons. 

In this review, we have presented a status report of our broad and systematic 
study of QCD-instanton induced processes in deep-inelastic scattering f| §, |10|, |ll 



n, 15 



We have emphasized that deep-inelastic scattering may be viewed as a distin- 
guished process for studying manifestations of QCD-instantons, since the available 
hard scale provides a dynamical infrared cutoff for the instanton size |12|]. At high 
Q 2 , instanton- induced fixed-angle processes may be reliably calculated, like in per- 
turbative QCD. Thus, deep-inelastic scattering at HERA offers a unique window to 
explore footprints of QCD-instantons. 

Of great interest is, of course, a first estimate of the total instanton-induced cross- 
section at HERA. We pursued a complementary approach to Ref. J7|, by working 
out ]14| a momentum-space picture of the instanton-contribution to the parton struc- 
ture functions and further inclusive observables. In leading-order of the semi-classical 
approximation about the instanton- antiinstanton configuration and at large Q 2 , the 
instanton-contribution to the gluon structure function J z 2g has the form of a convo- 
lution of a "splitting function" in the instanton-background with a total cross-section 
for the instanton-induced q*g subprocess. The latter contains the essential instanton 
dynamics. We presented a preliminary theoretical estimate of the total instanton- 
induced cross-section for deep-inelastic scattering at HERA, subject to appropriate 
kinematical cuts. Notably, the Bjorken scaling variable of the q*g subprocess has to be 
restricted from below, in order to retain theoretical control. Being in the (9(1 — 100) pb 
range, the resulting cross-section is surprisingly large. Due to inherent uncertainties 
associated with the renormalization- and factorization-scale dependencies, which are 



presently being investigated |1J|, this cross-section is, however, still to be considered 
preliminary. 

On the phenomenological side, we systematically explored the discovery potential 
for instanton-induced events at HERA, by studying the characteristics of the final 
state: A current-quark jet along with a densely populated hadronic "band" of width 
At] = ±0.9 in the (??i a b, 0iab)-plane, reflecting the isotropic decay of an instanton- 
induced fireball. Characteristic features include the large total transverse energy, 
Et = O(20) GeV, the large multiplicity, nb an d = 0(25), and the flavour-democratic 
production of hadrons, leading in particular to abundant production of K° mesons. At 
present, we consider the characterization of instanton-induced events by these final- 
state features to be more robust than predictions based on cross-section estimates. 

From a recent measurement of K° production in deep-inelastic scattering at 
HERA, the HI collaboration placed a first upper limit of 0.9 nb (95% confidence 
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level) on the cross-section for instanton-induced events in the kinematical region 
considered. This limit is less than an order of magnitude above our (preliminary) 
theoretical estimate, which even contains stronger additional cuts. We have studied 
in detail, how a combination of event shape information with searches of K° mesons, 
muons, and multiplicity cuts may further help to discriminate the QCD-instanton 
induced processes from the standard perturbative QCD background. 
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